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ture-sensing  probes  HUV.  HaTris,  NRL  Report  75B7,  Jtriy  -iW]  made  of  2^*m-diameter  platinum 
wire  which  were  separated  10  cm  and  optimally  excited  at  10 #WT  The  RMS  value  of  the  temper¬ 
ature  difference  between  the  probes  was  determined  with  an  R.  W.  Harris  double-AC-bndge/iog-metcr 
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20  ABSTRACT  (Continued) 

rand  the  meteorological  parameters  comprising  air  temperature,  dew  point,  wind  speed  and  direction, 
solar  insolation,  and  barometric  pressure  were  sampled  at  1-second  intervals.  Results  indicate  a  height 


dependence  of  the  .form  0jj^=  KH ",  where  K  is  a  function  of  the  surface  (0.84  <  K/l  O' 

1.8)  and  M  =  -1.2±0.5,  with  the  wide  variation  in  M  due  to  changing  meteorological  conditions,  most 
importantly  solar  insolation  and  wind  speed.  In  this  experiment,  no  significant  differences  in  aver¬ 
age  value  of  Pj^could  be  determined  for  the  lour  surfaces.  Additionally,  the  anticipated  reduction  in 
turbulence  over  aluminized  Mylar  was  not  observed 
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ATMOSPHERIC  TURBULENCE  AS  A  FUNCTION  OF 
HEIGHT  ABOVE  VARIOUS  SURFACE  MATERIALS 

INTRODUCTION 


This  report  presents  the  results  of  an  experiment  designed  to  determine  the  turbulence  strength 
over  four  surfaces  having  differing  reflective  and  physical  properties.  The  objective  is  to  use  these 
results  in  selecting  a  method  of  treating  surfaces  near  a  high-energy  laser  beam  which  yields  minimum 
thermally  induced  turbulence  and  reduced  optical  degradation  of  the  high-powered  beam 

Light  waves  propagating  through  the  atmosphere  exhibit  random  fluctuations  in  amplitude  and 
phase  due  to  variations  in  the  index  of  refraction  along  the  path.  An  example  of  this  is  starlight  scintil¬ 
lation.  Refractive-index  variations  are  directly  related  to  pressure,  temperature,  and  humidity  inhomo¬ 
geneities  resulting  from  turbulent  motion  in  a  real  atmosphere,  with  the  turbulent  motion  being  due  to 
changes  in  net  radiative  flux,  wind  speed  and  direction,  and  surface  roughness.  Much  theoretical  and 
experimental  work  has  been  performed  to  attempt  an  understanding  of  propagation  through  turbulent 
regimes.  Specifically,  the  efforts  of  Kolmogorov  [1),  Obukov  [2],  Chernov  13 ] ,  and  Tatarskii  [4,51 
should  be  mentioned;  much  of  the  theory  required  for  understanding  the  research  documented  in  this 
report  is  found  in  these  references. 


The  inhomogeneities  produced  by  atmospheric  turbulence  are  characterized  by  an  inner  tur¬ 
bulence  scale  /o,  representing  the  smallest  eddy  size  that  can  be  sustained  against  viscous  losses,  and  an 
outer  turbulence  scale  L0,  representing  large-scale  effects  resulting  from  wind  velocity,  wind  shear, 
height  above  the  ground  or  water,  and  the  like.  The  inner  scale  is  usually  about  1  to  10  mm,  depend¬ 
ing  on  turbulence  strength;  the  more  developed  the  turbulence,  the  smaller  the  inner  turbulence  scale. 
The  outer  turbulence  scale  can  be  about  1  to  10  m  and  represents  the  largest  scale  for  which  turbulence 
is  homogeneous  and  isotropic. 


The  random  variations  in  the  index  of  refraction  can  be  described  by  a  structure  function: 

D„(t)  *  (n(x)  -  nix  +  r)P.  (1) 

where  r  is  the  displacement  between  two  points  in  space.  A  similar  structure  function  can  be  defined 
for  temperature: 

Dt(  r)  =  [  TUr-n^JV .  (2) 

Under  the  assumption  of  locally  isotropic  and  homogeneous  turbulence,  Kolmogorov  derived  a  struc¬ 
ture  function  in  a  form  now  commonly  referred  to  as  the  two-thirds  law: 

Dn(r)  =  C„2  r2/}  (3) 

and 


DT(r)  «  Cj  r2li. 


(4) 


where  /n  <  r  =  |r|  <  L„  and  where  C„  and  C,  are  respectively  the  structure-function  parameters  for 
refractive  index  and  temperature.  From  Eqs.  (3)  and  (4)  one  can  see  that  the  turbulence  strength  is 
reflected  in  the  value  of  C'„2  or  Cj.  The  structure-function  parameters  are  related  by  the  expression  [6] 


C„  = 


11. bP 


1  + 


0.00753 


io-hc7 


(5) 
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where  P  is  the  pressure  in  millibars,  T  is  the  temperature  in  kelvins,  and  X  is  the  wavelength  in 
micrometers.  For  visible  wavelengths,  such  as  those  of  a  0.633-pm  He-Ne  laser,  Eq.  (5)  becomes 


C.  = 


19P 


10  bC,. 


(6) 


the  form  usually  seen  in  the  literature  [71.  By  combining  Eqs.  12)  and  (4)  and  substituting  the 
resulting  expression  for  C/  into  Eq.  (6),  one  arrives  at  an  equation  relating  variations  in  index  of 
refraction  to  a  temperature  differential: 


C„  = 


797 


io- 


(r,  -  r,)2 


1/2 


-1/3 


(7) 


where  Tt  and  T2  are  the  temperatures  at  two  points  separated  a  distance  r.  Equation  (7)  will  be  used  to 
determine  C„2  from  measurements  of  the  RMS  value  of  the  temperature  difference  A  T  between  two 
fine-wire  fast-response  thermal  probes,  as  has  been  done  previously  [6, 8, 9). 


The  structure  parameter  C„2  can  also  be  determined  from  optical  scintillation  measurements  [9); 
however,  this  measurement  requires  double-ended  experiments  over  paths  100  m  long  if  reasonable 
signal-to-noise  ratios  are  to  be  obtained.  For  the  present  work,  the  availability  of  a  mobile  thermal- 
probe  system  and  the  requirement  to  study  several  different  surface  materials  indicated  that  deriving  C„2 
from  measurements  of  A  T  would  be  more  feasible  in  a  limited  time  frame  and  much  less  expensive 
than  the  use  of  optical  scintillation  techniques. 


EXPERIMENT  AREAS  AND  EQUIPMENT 


Sites  and  Surfaces 


During  October  and  November  1977,  measurements  were  made  of  CT  over  open-field  grass, 
blacktop,  concrete,  and  aluminized  Mylar  at  the  Chesapeake  Bay  Detachment  (CBD)  of  NRL.  The 
open-field  grass  comprises  a  relatively  extensive  area  of  open  meadow  with  gradual  slopes  in  what  is 
called  West  Field.  Measurements  over  the  other  three  surfaces  were  performed  next  to  Building  49  in 
another  section  of  CBD,  somewhat  closer  to  buildings  and  trees.  This  second  experimental  site  con¬ 
tained  blacktop  and  concrete  areas.  In  addition,  32  sheets  of  4-by-8-ft  plywood,  covered  with  alumin¬ 
ized  Mylar  and  stapled  on  the  underside,  were  placed  on  the  concrete  area  to  make  a  95-m2  (1024-ft2) 
extended  squ,.,e  area,  as  sho>  n  in  Fig.  1.  Figure  2  provides  the  dimensions  of  the  blacktop  and  Co.,- 
crele  areas  and  shows  the  lot..' inn  of  the  trailer  van  which  housed  the  data-acquisition  system  and  the 
electronics  for  the  meteorological  sensors  and  thermal  probes. 

Mobile  Tower,  Meteorological  Sensors,  and  C'j  Probes 

Figures  3  and  4  show  the  mobile  tower  on  which  sensors  were  mounted  for  monitoring  air  tem¬ 
perature,  dew  point,  wind  speed,  and  wind  direction.  Also  shown  in  Fig.  4  are  the  fast-response  ther¬ 
mal  probes  used  for  measuring  C;  and  the  radiometer  used  for  determining  global  radiation. 


Air  temperature  and  dew  point  were  monitored  with  an  ECi&G  model  ilOS-M  automatic  dew¬ 
point  hygrometer  having  a  range  of -62°C  to  49°C',  a  response  time  of  1.7°C/s,  and  an  accuracy  of 
0  3°C  Wind  speed  and  direction  were  obtained  from  an  R.  M.  Young  Company  model  35003  Gill  pro¬ 
peller  vane  having  a  wind-speed  threshold  of  0.1  to  0.2  m/s  and  a  range  of  0  to  30  m/s.  Global  radia¬ 
tion.  the  sum  of  direct  solar  radiation  and  diffuse  atmospheric  radiation,  was  measured  using  an  Eppley 
Laboratory  thermopile  pyranometer  having  a  response  time  of  3  to  4  s  and  an  output  of  1 1.25  V  per 
W/m2.  Barometric  pressure  was  monitored  with  a  Yellow  Springs  Instrument  Company  Model  2014 
pressure  transducer  having  a  range  of  833  to  1067  millibars;  this  sensor  was  located  inside  the  trailer 
van. 
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l  \  -  Mobile  lower  with  meteorological  sensors  deployed  above 
aluminized-Mylar  covered  plywood  sheets 


|  ,s  2  _  Layout  of  concrete  and  blacktop  surfaces  available 
for  turbulence  measurements 


Fig.  3  —  Mobile  tower  with  sensors  for  monitoring 
wind  direction,  wind  speed,  Cr,  air  temperature, 
and  dew  point  (front  view) 
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Figure  5  shows  the  upper  section  of  the  tower.  At  the  top  is  the  Young  wind  speed  and  direction 
sensor;  on  the  right-hand  side  of  the  tower  is  the  EG&G  sensor.  The  left-hand  side  of  the  tower  shows 
the  welded  section  of  Unistrut  which  supported  the  movable  aluminum  L  section.  The  fast-response 
thermal  probes  were  attached  to  this  section,  which  was  bolted  in  place  at  the  desired  height.  For 
heights  over  3  m  the  entire  tower  was  tilted  to  a  horizontal  position  to  facilitate  moving  the  C,  probes 
to  a  different  level. 


The  dual  fine-wire,  fast-response  temperature  probes  used  for  determining  C„2  through  the  meas¬ 
urement  of  Cr  are  shown  in  Fig.  6.  The  vane  assembly  kept  the  two  probes  pointed  into  the  wind. 
F.ach  probe  uses  platinum  wire  2.032  gm  in  diameter  and  1.0  mm  in  length,  with  the  separation 
between  the  two  probes  being  10  cm.  These  probes  are  of  the  same  design  as  that  used  by  Dowling 
and  Livingston  19].  They  are  part  of  an  AC  bridge  which  allows  measurement  of  the  temperature  of 
each  probe.  The  probes  are  optimally  excited  with  10giW  of  electrical  power  in  accordance  with  the 
results  of  Harris  [10],  The  electronics  required  for  the  thermal  probes  are  divided  into  a  front-end  sec¬ 
tion  positioned  near  the  probe  assembly  (at  the  base  of  the  tower)  and  a  back-end  section,  containing  a 
B&K  Instruments  model  215  RMS-voltmeter/log-converter,  installed  in  the  trailer  van.  The  entire  sys¬ 
tem  was  originally  designed  by  R.  W.  Harris;  a  redesign  of  the  front-end  unit  and  modifications  of  the 
back-end  unit  have  been  performed  by  Tower  Systems,  Inc.  Similar  fast-response,  thermal-probe  sys¬ 
tems  are  described  in  Ref.  8.  Figure  7  shows  one  of  the  fine-wire  probes  with  the  protective  cap  on  and 
the  other  with  the  cap  removed. 

All  of  the  voltages  from  the  previously  described  sensors  were  processed  by  a  Monitor  Labs  7200 
data  logger  which  was  connected  to  a  Digidata  130  Series  nine-track,  incremental,  digital  tape  recorder. 
The  digital  tapes  were  processed  on  the  Texas  Instruments  Advanced  Scientific  Computer  at  NRL  using 
Fortran  computer  programs  developed  for  this  purpose.  Computer  plots  of  the  data,  presented  in  a 
later  section,  were  also  obtained.  Figure  8  is  a  diagram  of  the  sensors  and  data-acquisition  system. 
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Dual  fine-wire,  fast -response  temperature  probes  lor  measuring  ( 
(shown  mounted  on  a  wind -vane  assembly) 


7  —  line-wire  thermal  probe  with  the  protective  cap  (top) 
and  a  probe  with  the  protective  cap  removed  (bottom) 
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EXPERIMENTAL  PROCEDURE 

Turbulence  levels  were  measured  over  grass  in  West  Field  and  over  blacktop,  concrete,  and 
aluminized  Mylar  next  to  Building  49  at  CBI).  The  experimental  procedure  followed  was  identical  for 
each  of  the  four  surfaces.  The  mobile  tower  was  placed  such  that  the  end  where  the  turbulence  probes 
were  lot;  ed  was  pointed  into  the  prevailing  wind  direction  to  maximize  the  effect  of  the  surface  being 
studied  and  minimize  the  contribution  to  the  turbulence  measurement  resulting  from  mechanical  mix¬ 
ing  of  the  air  surrounding  the  tower  structures.  The  probes  were  installed  on  the  vane  assembly,  the 
protective  caps  on  each  probe  were  removed,  and  the  protective  sheath  on  each  probe  was  slid  back  to 
its  farthest  position.  The  exposed  probes  were  left  in  position  overnight,  if  needed,  when  weather  con¬ 
ditions  indicated  no  near-term  precipitation,  since  the  platinum-wire  structure  is  relatively  sturdy  even 
in  high  winds.  The  Fpplcy  radiometer  was  taken  out  of  the  trailer  van  each  day  and  placed  outside  next 
to  the  tower.  The  turbulence  front-end  unit,  requiring  117-V  AC  power,  was  installed  on  the  tower 
trailer  each  day,  since  this  unit  was  not  weatherproof. 

All  of  the  meteorological  sensors  and  the  data-acquisition  systems  were  activated  1/2  hour  prior  to 
the  taking  of  data  The  time  constant  for  averaging  the  temperature  difference  between  the  two  probes 
was  set  to  1  s  on  the  B&K  RMS-voltmeter/log-converter.  Thus  the  averaging  period  of  Eq.  (6)  is  the 
time  constant  of  the  B&K  instrument  The  Monitor  Labs  data  logger  scanned  all  input  channels  once 
every  second. 

Two  sizes  of  computer  tape  were  used  a  183-m  (600-ft)  reel  and  a  366-m  (1200-ft)  reel,  depend¬ 
ing  on  the  anticipated  length  of  a  particular  day's  run.  Each  data  tape  was  processed  at  NRL  at  the  end 
of  the  day  or  on  the  next  day  to  insure  validity  of  the  data  and  to  check  on  the  performance  of  the  sen¬ 
sors  and  the  data-acquisition  system 

Measurements  of  (  t  were  made  at  heights  of  0.5,  1,  2,  3,  4.  and  5  m  above  the  surface.  For  the 
positions  from  0.5  to  3  m  on  the  tower,  the  turbulence  probe  was  moved  to  the  level  and  bolted  in 
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place;  the  4-  and  5-m  positions  required  that  the  tower  be  tilted  to  the  horizontal  before  moving  the 
probe  assembly  to  the  next  position.  This  latter  procedure  took  approximately  5  to  10  min.  Data  were 
taken  at  each  level  for  intervals  ranging  from  15  to  45  min.  Each  minute  yielded  60  data  points  of  the 
meteorological  parameters  and  of  C it  was  expected  that  the  approximate  1/2-hour  data  period  at  each 
level  would  suffice  for  determining  ( for  that  height  above  a  particular  surface. 

Data  were  usually  taken  during  the  hours  0900  to  1500.  Changes  in  the  weather  required  shorten¬ 
ing  the  interval  during  which  data  were  taken,  especially  when  rain  was  expected  or  occurred.  One  run 
was  made  from  1630  on  14  November  until  0730  the  next  day  to  determine  the  overnight  variation  of 
C;  at  a  level  of  3  m  over  an  aluminized-Mylar  surface. 

For  measurements  over  blacktop  (Fig.  2)  the  mobile  tower  was  placed  at  the  northwest  corner  of 
the  concrete  area  with  the  probes  pointed  into  a  northwesterly  wind.  The  measurements  over  the  rec¬ 
tangular  concrete  area  (Fig.  2)  were  made  with  the  mobile  tower  at  the  center  with  the  probes  pointed 
into  a  westerly  wind. 

The  32  plywood  sheets  covered  with  aluminized  Mylar  were  placed  in  front  of  the  mobile  tower 
each  morning  and  then  taken  up,  covered,  and  stored  in  Building  49  each  evening  (except  for  the  one 
overnight  run).  Through  this  procedure  the  mirrorlike  reflectivity  of  this  material  was  maintained  dur¬ 
ing  the  experimental  period. 

RF.SIT.TS 

Cn  Over  Crass 

Figures  9,  Hi.  I  1,  and  13  arc  each  a  computer  plot  of  the  refractive-index  structure  parameter  C,; 
-denoted  (  NSQ).  solar  radiation,  and  wind  speed  versus  time  of  day  derived  from  measurements  over 
cr  iv-  I  he  refractive- index  structure  parameter  is  plotted  on  a  logarithmic  scale  as  a  solid  line.  Solar 
..’.tuition,  which  includes  all  sol., r- derived  radiation  (previously  defined  as  global  radiation),  is  plotted 
in  a  linear  scale  as  a  dotted  line.  Wind  speed  is  similarly  plotted  on  a  linear  scale  as  a  dashed  line  The 
-’at. i  points  plotted  arc  3-mm  averages,  which  in  most  instances  represent  the  average  of  ISO  measure¬ 
ments  obtained  using  the  l-s  scan  period  of  the  data  logger 

The  measurements  made  on  12  October  1 977.  are  shown  in  Fig.  9.  The  skies  were  overcast,  and 
the  wind  was  out  of  the  northwest  at  about  2  to  5  m/s.  as  shown  by  the  dashed  curve.  The  dotted 
solar-radiation  line  showed  a  general  increase  during  the  day;  the  variations  are  manifestations  of  the 
changing  cloud  cover. 

The  two  vertical  bars  above  the  axis  show  the  last  data  points  for  the  1.3-m  and  3.7-m  heights  at 
which  the  („'  measurements  were  made.  The  value  of  shown  by  the  solid  line,  was  influenced 
continuously  by  the  nature  of  the  surface,  by  the  solar  radiation,  and  by  the  wind  speed  in  a  complex 
interactive  process.  It  is  this  aspect  which  makes  analysis  difficult  and  comparisons  between  days  sub¬ 
ject  to  varying  interpretation.  The  following  comments  on  the  data  therefore  will  tend  to  be  more 
qualitative  than  quantitative. 

Figure  9  shows  three  important  effects:  the  influence  of  wind  speed  on  C„\  the  close  relationship 
between  (  and  solar  radiation,  and  the  decrease  in  turbulence  at  greater  heights.  The  first  of  these 
effects  is  seen  at  1227.  From  1224  to  1230  (he  solar  radiation  was  approximately  constant  at  0.65 
W/nr,  and  the  height  H  above  the  grass  was  3.7  m  Thus  two  of  the  important  influences  on  C„2  were 
essentially  constant  The  distinct  peak  in  wind  speed  produced  an  equally  distinct  trough  in  C„\  leading 
to  the  conclusion  that  increased  wind  speed  causes  decreased  small-scale  turbulence  through  the  injec¬ 
tion  of  relatively  homogeneous  masses  of  air  into  the  measurement  region. 
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F-ig.  9  —  ( „3  nbovc  grass  (solid  line!,  solar  raduuion  (dolled  line! 

and  wind  speed  (dashed  line)  on  12  October  1977 


The  seeond  important  effect  is  seen  throughout  the  day:  peaks  in  the  C„2  curve  correlate  with 
peaks  in  the  solar  radiation  curve,  and  troughs  in  ( correlate  with  troughs  in  solar  radiation.  Exam¬ 
ples  of  this  occur  at  1042  (peaks  in  Or  and  solar  radiation)  and  at  1051,  1106.  1121  and  1209 
(troughs).  If  one  notes  that  from  10.10  to  114'  the  solar  radiation  has  an  average  value  of  0.4  W/nr, 
one  can  estimate  the  influence  of  solar  radiation  on  the  structure  parameter,  the  drop  to  0.27  W/nr  in 
solar  radiation  induces  a  drop  ip  from  its  nominal  value  of  X  *  10  |J  to  5  x  10  11  m  :  \  a  drop  of 
almost  50%.  In  conjunction  with  this  correlation  of  solar  radiation  and  one  should  note  that 
changes  in  the  structure  constant  occur  with  a  time  lag  of  about  2  to  6  min.  in  contrast  to  the  time  lag 
for  the  influence  of  the  wind  speed,  which  seemed  to  be  imperceptible  w  ith  the  3-min  resolution  of  the 
data  points. 

Finally,  the  effect  of  moving  the  probe  to  a  different  height  can  be  observed:  At  1202  the  probe 
assembly  was  moved  from  a  height  H  •“  1.3  m  (4  ft)  to  //  =  3.7  m  (12  ft),  a  factor-of-three  increase. 
Between  1209  and  1300  the  value  of  („2  has  a  slightly  lower  average  of  7  x  10  14  m  2  ’.  even  though 
the  solar  radiation  has  a  net  increase  of  0.2  W/m2  per  hour.  This  observation  implies  that  the  effect  of 
solar  radiation  on  C’,,2  is  lessened  at  heights  above  3  m  and,  most  importantly,  that  tripling  the  height 
above  a  grass  surface  essentially  counteracts  the  increase  in  („2  that  would  occur  with  the  increase  of 
solar  radiation  from  an  average  0.4  W/m2  to  0.6  W/m2  between  the  heighis  13  and  3.7  m.  a  50% 
increase  in  solar  radiation.  Had  the  solar  radiation  remained  constant  during  the  measurements  at  these 
two  heights,  one  would  expect  that  („2  would  be  lower.  This  agrees  with  Wyngaard  et  al.  |7),  whose 
own  data  agree  with  an  H  4/J  dependence  of  C’j  for  local  free  convection  on  a  summer  day  with  light 
winds  at  heights  a  few  meters  above  the  ground,  an  experimental  situation  similar  to  the  one  of  this 
report,  except  for  the  overcast  conditions.  To  eliminate  the  effects  of  changing  wind  speed  and  solar 
radiation,  multiple  probe  assemblies  at  several  heights  would  be  required  to  deduce  the  functional 
dependence  of  ( „2  on  //.  (The  use  of  the  one  available  probe  system  made  it  more  difficult  to  deter¬ 
mine  an  experimental  value  for  the  exponent  of  H.) 

Figure  10  shows  the  results  of  the  measurements  on  17  October  1977.  at  heights  of  0  5.  10,  and 
2.0  m.  The  day  was  one  of  high  winds,  from  7  to  14  m/s,  and  clear  skies.  The  inverse  relationship 
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tig  10  —  <  „!  above  grass  (solid  line*,  solar  radiation  (dolled  line). 

and  wind  speed  (dashed  line)  on  17  October  1977 

between  C„2  and  wind  speed  is  again  evident  at  1221,  1242,  1336,  1412,  1418,  1424,  and  1447,  though 
an  increase  in  C„2  to  3  x  10" 13  m  2/3  at  1351  seems  to  have  been  caused  by  an  increase  in  wind  speed  to 
a  very  high  13.7  m/s,  an  effect  opposite  to  what  had  been  usually  observed. 

From  1215  until  1400  the  solar  radiation  dropped  very  little  and  was  essentially  0.75  W/m2;  dur¬ 
ing  this  time  the  winds  were  gusty,  with  an  average  value  of  10  m/s.  For  H  =  0.5  m,  C„2  -  4.5  x 
10  13  m"223;  for  H  =  1.0  m,  C'„2  -  1.9  x  10"13  m~v\  for  H  =  2.0  m,  C„2  =  0.9  x  10  13  m  2/J.  By  plot¬ 
ting  these  three  data  points  on  a  log-log  format,  one  can  obtain  the  relationship  between  C„2  and  height 
H  above  the  grass  on  this  day:  C'„2  —  1.9  x  I0~13  H  1  IAt008  m"2/3.  The  theory  for  unstable  free  con¬ 
vection  with  no  wind  shear  indicates  an  W~4/3  dependence,  stable  conditions  indicate  an  H  Vi  depen¬ 
dence.  Thus  the  experimentally  determined  exponent  -1. 16+0.08  falls  in  the  proper  range  between 
-2/3  and  -4/3;  that  it  is  closer  to  -4/3  satisfies  the  initial  assumption  of  the  occurrence  of  free  con¬ 
vection  (7],  The  drop  in  solar  radiation  after  1400,  which  corresponds  to  the  data  for  2.0  m,  would 
tend  to  decrease  the  values  of  C„2  during  this  time;  this  tends  to  increase  the  negative  slope  on  the  log- 
log  plot,  thus  indicating  that  either  the  exponent  lies  in  a  region  close  to  —1  or  no  simple  exponential 
relationship  of  the  form  derived  in  the  theory  can  describe  what  occurs  in  the  region  below  2  m. 

Figure  1 1  shows  the  results  of  the  measurements  on  18  October  1977.  No  solar-radiation  data  are 
plotted,  because  of  a  malfunction  that  day.  Since  the  skies  were  clear,  one  would  expect  an  increase  in 
solar  radiation  from  900  to  1200  and  a  gradual  decrease  thereafter.  The  winds  were  out  of  the 
southwest  at  about  3±  1.5  m/s.  Figure  12  is  a  log-log  plot  of  C'„2  versus  height  for  the  data  shown  in 
Fig.  1 1.  The  line  drawn  has  the  form  (  „2  =  1.8  x  10  13  H  1  72  m  2/3,  determined  from  a  graphical  esti¬ 
mate  of  where  the  line  should  be.  The  power  of  //  is  greater  than  the  -4/3  predicted  by  free- 
convection  turbulence  theory.  The  influence  of  solar  radiation  is  most  apparent  for  the  0.5-rn  data  (Fig. 
12).  for  which  decreased  sunlight  produced  a  value  of  C„2  lying  below  the  line. 

Figure  13  shows  the  measurement  results  for  19  October  1977,  at  heights  of  0.5,  1.0,  and  2.0  m. 
The  respective  average  values  of  (  „2  at  these  heights  are  1.5  x  10  l3,  0.8  x  10  l3,  and  0.5  x  10  13 
m  2  these  yield  the  relation  (  =  8.4  x  10  14  //  11  sl  m  2'  \  which  has  a  power  of  //  quite  different 
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Fig  1 1  —  C'rtJ  above  grass  (solid  line),  and  wind  speed 
(dashed  line)  on  18  October  1977 


Pig  12  —  Dependence  of  (  „*  (above  grass)  on  height 
for  18  October  1977 
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from  that  obtained  on  the  previous  day.  However,  the  rapid  increase  in  solar  radiation  during  the 
measurement  interval  (0.22  W/nv  per  hour)  results  in  large  C,2  values  for  1.0  and  2.0  m  when  com¬ 
pared  with  the  case  of  constant  solar  radiation  This  changing  solar  radiation  does  not  allow  an  unambi¬ 
guous  determination  of  the  variation  of  C  with  height  for  these  data 

t'n  Over  Blacktop 

Measurements  of  C,j  were  made  over  blacktop  on  20  October  1977.  The  results  are  presented  in 
Fig.  14.  The  three  effects  noted  over  grass  can  be  seen  again:  First,  at  1012  and  1218  the  dip  in  wind 
speed  produces  a  concomitant  peak  in  („2,  even  though  these  times  correspond  to  the  disparate  5.5-  and 
1.0-m  heights.  Second,  the  large  drop  in  solar  radiation  occurring  at  1306  is  paralleled  by  a  significant 
drop  in  C'„2.  Finally,  as  the  measurement  height  above  the  blacktop  decreased,  the  value  of  the  struc¬ 
ture  parameter  increased.  For  1.0  ,  2.0-.  and  4.0-m  heights,  the  respective  average  values  of  C'„2  are  1.2 
x  10" '  \  4.2  x  10  M,  and  2  0  x  10  14  m  2  \  a  log-log  plot  of  these  data  yields  C„2  =  1.2  x  10  "  H  ' 32 
m  2' \  which  generally  has  the  same  constant  and  power  of  H  as  previously  found  in  measurements 
over  grass.  It  is  surprising,  however,  to  find  that  //  has  this  nearly  -4/3  exponent  even  though  the 
solar  radiation  increased  throughout  most  of  the  measurement  period.  If  the  same  reasoning  is  used 
here  as  was  used  for  the  October  19  grass  measurements,  compensation  for  the  effect  of  increasing 
solar  radiation  would  lessen  the  values  of  C„2  at  2.0  and  1.0  relative  to  that  at  4.0,  thus  decreasing  the 
power  of  H,  that  is,  yielding  a  value  closer  to  1.0.  This  implies  that  turbulence  above  blacktop  tends  to 
decrease  less  rapidly  with  height  than  it  does  over  grass,  which  might  be  expected  for  a  surface  which 
absorbs  more  solar  energy  and  becomes  hotter. 

C„2  Over  Concrete 

Measurements  were  made  on  31  October  1977,  over  concrete,  and  the  results  are  shown  in  Fig. 
15.  Approximate  average  values  of  C„2  at  0.5,  1.0,  2.0,  3.0,  4.0,  and  5.0  m  are  respectively  3.2,  2.0, 
4.5,  3.0,  2.3,  and  2.0  x  I0"14  m2/\  The  day  became  cloudy  at  around  1330,  at  which  time  solar  radia¬ 
tion  and  wind  speed  fell  rapidly;  measurements  at  1.0  m  were  ongoing  during  this  change  in  conditions. 
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The  C2  values  for  0.5  and  1.0  m  are  approximately  an  order  of  magnitude  lower  than  they  would  have 
been  had  the  solar  radiation  remained  at  0.75±0.05  W/m2,  that  for  2.0  to  5.0  m.  The  changes  in 
weather  conditions  essentially  change  the  multiplicative  constant  in  the  C 2  equation.  The  data  lead  to 
two  equations: 

C„2=  2.0  x  I0‘u/r068  m~2/3.  H  =  0.5,  1.0  m; 

and 

C2  =  9.0  x  10  14  H -  °m ~2J\  H  =  2.0  to  5.0  m. 

These  two  equations  are  derived  from  Fig.  15;  the  values  of  C2  at  the  lower  heights  are  small,  thus 
leading  to  large  differences  in  the  structure  constant  at  these  levels,  with  a  resulting  greater  uncertainty 
in  the  exponent  of  H.  Figure  15  at  first  would  seem  to  indicate  lower  C2  values  for  concrete  than  for 
blacktop  (Fig.  14).  However,  for  the  blacktop  data  the  solar  radiation  was  increasing  during  most  of  the 
measurements,  whereas  for  the  concrete  the  solar  radiation  was  essentially  unchanged  for  H  above  3.0 
m.  The  two  previous  equations  for  C2  actually  give  values  for  the  structured  constant  which  are  larger 
than  those  calculated  from  the  blacktop  C„2  equation.  The  intuitive  concept  that  turbulence  levels 
above  concrete,  which  tends  to  reflect  more  solar  radiation,  should  be  less  than  those  over  blacktop 
thus  is  not  verified.  Perhaps  this  is  due  to  the  uncertainties  in  the  data  for  concrete;  a  more  likely 
explanation  is  that  the  experimental  concrete  area  was  rather  small  and  also  surrounded  on  three  sides 
by  blacktop,  with  turbulence  over  the  blacktop  influencing  the  turbulence  over  the  concrete. 

Over  Aluminized  Mylar 

Figures  16  through  19  show  the  results  of  measurements  made  over  an  aluminized-Mylar  surface. 
The  data  for  1  November  1977,  shown  in  Fig.  16  were  obtained  using  only  12  of  the  32  sheets  of 
aluminized-Mylar-covered  plywood.  The  solar  radiation  increased  from  0951  to  1115  and  decreased 
thereafter  due  to  the  increasing  cloud  cover.  The  effect  of  broken  cloud  cover  prior  to  1130  is  evident 
The  winds  remained  at  about  3  m/s.  The  variations  in  solar  radiation  make  it  difficult  to  extract  quanti¬ 
tative  changes  in  C'2  at  the  various  heights,  though  the  drop  in  turbulence  with  increasing  height  is 
apparent. 


:  i  »-*i 

J  ig  16  —  (  *  above  alumini/ed  Mvlar  (solid  line),  solar  radiation  (dotted  line), 
and  wind  speed  (dashed  line)  on  l  November  1977 
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Figure  17  shows  the  results  for  11  November  1977,  a  day  beset  by  broken  clouds  with  winds  2±  1 
m/s  out  of  the  west  to  southwest.  The  broken  cloud  cover  produced  wide  variations  in  solar  radiation 
which  led  to  large  changes  in  f,2;  the  effect  is  clearly  evident  at  1  103,  1148,  1212,  1233,  and  1336 
(peaks)  and  at  1  142,  1203,  1318,  and  1354  (troughs).  Because  of  these  many  variations  in  solar  radia¬ 
tion,  it  is  again  difficult  to  determine  quantitatively  the  height  dependence  of  C„2;  however,  the  general 
tendency  to  increased  C~  with  decreased  height  is  again  apparent.  More  importantly,  this  figure  shows 
that  aluminized  Mylar  produces  turbulence  levels  which  have  a  wider  response  to  solar  radiation,  rang¬ 
ing  from  5.25  x  10"  m  2/3  at  1 142  to  3.44  x  10"  m  2/3  at  1412. 


I  it  I'  -  <  .tht'VL-  .ilumini/cil  Mvl.u  Nnlul  1 1 no > .  xnl.ir  rjdialton  Uloltod  lino), 
jnd  wind  speed  Mashed  line*  on  1  I  November  W7 


f  igure  18  shows  the  results  of  data  taken  on  14  November  1977.  Some  clouds  caused  variations 
in  the  solar  radiation  during  the  measurement  period.  The  direct  relation  between  C„2  and  solar  radia¬ 
tion  is  evident  in  the  troughs  at  1009,  1036,  1124,  1200,  and  1218  and  in  the  peaks  at  1151,  1209,  and 
1224.  Estimated  average  C„2  values  at  0.5.  1.0,  and  2.0  m  are  respectively  2.8,  1.4  and  0.86  x  10'  13 
m2  \  which  lead  to  (  *  1.5  x  10  13  //  m  2  3  .  The  influence  in  the  power  of  H  due  to  increasing 

solar  radiation  is  evident;  had  the  solar  radiation  remained  constant,  the  C„2  value  at  1.0  and  2.0  m 
would  have  been  lower,  producing  a  higher  value  for  the  exponent  of  //,  which  is  in  the  correct  direc¬ 
tion  and  perhaps  nearer  the  1  to  1.3  obtained  for  the  other  surfaces.  The  larger  variation  in  C „2  with 
change  in  solar  radiation  is  clearly  seen;  that  is,  the  turbulence  level  above  aluminized  Mylar  is  more 
sensitive  to  solar  radiation  than  it  is  above  grass,  concrete,  and  blacktop. 

The  results  for  18  November  1977,  are  shown  in  Fig.  19.  The  skies  were  clear  until  1100,  when 
some  clouds  appeared,  causing  significant  drops  in  solar  radiation.  The  dependence  of  ( „2  on  solar  radi¬ 
ation  is  especially  evident  at  1103  and  1203.  Values  for  (  „2  at  1.0,  2.0,  3.0,  4.0,  and  5.0  m  are  respec¬ 
tively  18,  7.5,  5.0,  2.8,  and  2.3  x  10  14  m  2/\  which  lead  to  C„2  =  1.8  x  10' 13  H  1 3  m'2'3.  This  power 
of  //  is  in  agreement  with  the  theoretical  value  of  -4/3,  but  the  increase  of  solar  radiation  during  the 
measurements  would  imply  a  reduction  in  the  H  exponent. 
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For  the  four  measurement  days  with  aluminized  Mylar,  it  seems  that  the  height  dependence  of  („2 
is  proportional  to  a  power  of  H  in  the  neighborhood  of  —  1.  The  influence  of  solar  radiation  seems  to 
be  greater  than  for  the  prior  three  surfaces.  Originally  it  had  been  proposed  that  the  specular  nature  of 
aluminized  Mylar  would  lead  to  a  lowering  of  turbulence  above  this  surface  The  results  of  this  experi¬ 
ment  do  not  support  this  hypothesis;  they  showed  an  increase  in  the  dynamic  range  of  the  turbulence 
Again  turbulence  from  bordering  surfaces  may  be  significant. 

Overnight  Measurements  of  C„2  Over  Aluminized  Mylar  at  3.0  m. 

Figure  20  shows  the  results  of  measurements  made  from  1418  on  14  November  until  0730  on  15 
November  1977,  at  3.0  m  above  aluminized  Mylar.  The  skies  were  clear  during  the  night,  and  the 
winds  decreased  toward  sunset  and  rose  slightly  during  the  night.  The  logarithmic  scale  for  C„2  is 
compressed  by  an  additional  cycle  from  that  of  the  previous  figures  due  to  C„2  minima  below  10  15 
m~2/3.  Important  effects  that  can  be  seen  are  the  decrease  in  C„2  with  decrease  in  solar  radiation,  lead¬ 
ing  to  a  minimum  after  sunset,  when  the  air  temperature  and  surface  temperature  became  equal,  an 
increase  in  C„2  during  the  night,  due  to  the  reversal  in  direction  of  the  net  heat  flux  after  sunset  and  a 
mimimum  after  sunrise  when  the  air  temperature  and  surface  temperature  again  become  equal.  The 
increased  wind  speed  from  0300  to  0500  appears  to  reduce  the  turbulence  during  this  time,  in  accord 
with  the  results  obtained  for  the  daytime  measurements. 
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tig  20  —  (  „2  above  aluminized  Mylar  (solid  line),  solar  radiation  (dolled  line), 
and  wind  speed  (dashed  line)  on  14  and  IS  November  1977 


CONCLUSIONS 

Based  on  the  data  shown  in  Figs.  9  through  20,  the  following  observations  can  be  made: 

•  The  refractive-index  structure  parameter  C„2  is  inversely  related  to  wind  speed  at  heights 
up  to  5  m  over  the  four  surfaces  studied. 

•  C„2  is  a  direct  function  of  the  solar  radiation.  There  is  a  stronger  correlation  for  alumin¬ 
ized  Mylar  than  for  grass,  blacktop,  or  concrete. 

•  („2  is  an  inverse  function  of  the  height  above  the  surface;  the  relationship  has  the  form 

=  KH  w,  where  K  depends  on  the  type  of  surface,  0.84  4  A/10  13  m~2/3  <  1.8,  and 
M  varies  with  surface  and  meteorological  conditions.  For  this  experiment  0.68^  M  ^ 
1.72,  with  M  -  1.2  being  a  good  estimate  in  most  instances. 
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•  No  significant  difference  in  the  average  value  of  C„J  could  be  determined  in  this  experi¬ 
ment  for  the  four  surfaces  studied,  though  the  range  of  the  structure  parameter  was 
greater  for  aluminized  Mylar.  An  anticipated  reduction  in  Cj  over  the  aluminized-Mylar 
surface  was  not  observed  in  this  experiment. 

RECOMMENDATIONS 

Additional  insight  into  the  turbulence  processes  occurring  near  the  ground  could  be  obtained  by 
additional  measurements  if  two  modifications  were  made  to  the  present  experiment:  First,  expansion  of 
the  Cr  measurement  system  from  one  pair  of  fast-response  thermal  probes  to  perhaps  five  pairs  spaced 
I  m  apart;  this  would  allow  a  determination  of  the  height  profile  independent  of  wind  and  solar- 
radiation  variations.  Second,  use  of  larger  surfaces  would  ensure  that  the  differences  measured  between 
surfaces  are  due  to  the  surface  under  consideration  and  not  to  different  bordering  surfaces  and  struc¬ 
tures. 


The  results  presented  indicate  a  need  to  repeat  the  experiment  at  times  of  the  year  having  a 
greater  percentage  of  clear  days  with  higher  solar  insolation  and  less  wind,  such  as  occur  during  the 
summer. 
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